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Abstract 

In this article we propose a model for stochastic delay differential equation 
with jumps (SDDEJ) in a differentiable manifold M endowed with a connection 
V. In our model, the continuous part is driven by vector fields with a fixed delay 
and the jumps are assumed to come from a distinct source of (cadlag) noise, 
without delay. The jumps occur along adopted differentiable curves with some 
dynamical relevance (with fictitious time) which allow to take parallel transport 
along them. Using a geometrical approach, in the last section, we show that the 
horizontal lift of the solution of an SDDEJ is again a solution of an SDDEJ in 
the linear frame bundle BM with respect to a connection in BM. 
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1 Introduction 

Many natural phenomena present delays with respect to inputs: mathematically, 
this is a well known and established theory in the literature. In fact, standard delay 
differential equations have been extensively studied, see e.g. the classical Hale m 
and references therein. Models for these equations in manifolds appear in Oliva 
m, and stochastic perturbations are considered in Langevin, Oliva and Oliveira 
m, Mohammed [18| . Mohammed and Scheutzow EI1E2], Caraballo, Kloeden and 
Real [2] for SPDE, among many others. More recently, differential equations with 
random unbounded delay have been considered in Garrido-Atienza, Ogrowsky and 
Schmalfuss m- Besides delay, another usual characteristic in systems in biology, 
physics, economics, climatology, etc is the presence of jumps, both in the input and 
in the output. In this paper, we put these two characteristics, delays and jumps, 
together in the same mathematical framework. 

In our model of stochastic delay differential equation with jumps (SDDEJ), the 
continuous part of the solution is driven by vector fields with a fixed delay d > 0, 
and the jumps are assumed to come from a distinct source of (cadlag) noise, without 
delay. This idea is inspired by the fact that in several phenomena, informations reach 
a receptor by different communication sources (or channels), hence it is reasonable 
that delays on time are dependent on these distinct sources. As a simple example: 
in a storm, lightnings have instantaneous impact, but thunders come with delay. 

Following the ideas behind the so called Marcus equation, as in Kurtz, Pardoux 
and Protter |15] . our jumps in the solution occur along fictitious differentiable curves 
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which allow to take parallel transport along them. Generically speaking, as in Section 
2 , these fictitious differentiable curves can be taken in many distinct ways (say, 
randomly or along geodesics, etc). In our model, presented in Section 3, they follow 
the deterministic flow generated by the vector fields, without delay. We consider 
that the noise splits into a continuous component which is a Brownian motion with 
drift (extensible to a class of continuous semimartingales) plus a component given 
by a sequence of (cadlag) jumps. The number of jumps is assumed to be finite in a 
compact time interval. Hence it includes Levy-jump diffusion, but not Levy process 
in general. This idea of finite jumps in bounded intervals has parallel in the theory 
of chain control sets, see Colonius and Kliemann [5], Patrao and San Martin [20], 
and references therein. 

Another model with some numerical results for delay stochastic systems with 
jumps in Euclidean spaces can be found in Dareiotis, Kumar and Sabanis [7] for 
Levy processes. Also, stochastic geometry with jumps is considered in Cohen [3|, 
where the authors use second order calculus. 

In our approach, the delay are treated using parallel transport along the solu¬ 
tions, prescribed by a connection V in a differentiable manifold M. In Catuogno 
and RufRno |S|, the authors consider a geometrical approach to stochastic delay 
differential equations (SDDE) on a manifold M. In particular, they prove that the 
horizontal lift of a SDDE solution to the linear frame bundle BM is, again, a solution 
of an associated SDDE in BM, with respect to a horizontal connection in BM. 

The paper is organized as follows: In Section 2 we construct parallel transports 
along curves with jumps: here, the jumps are taken along a generic family of 
fictitious curves. In Section 3, our model of delay differential equations with jumps 
is presented. Compactness of the manifold M is assumed only to guarantee the 
existence and uniqueness of solution for SDDE (without jumps) for all t > 0, as in 
Leandre and Mohammed m- Finally, in Section 4, we explore geometrical aspects 
of the SDDEJ. After a short revision on the geometry of linear frame bundle and on 
the horizontal connection in BM, we show that the horizontal lift of solutions of 
SDDEJ in M are again solutions of SDDEJ in BM, extending the result mentioned 
above in |3|. 

2 General aspects of parallel transport 

Let M be a differentiable manifold, and V a connection on M. This structure, 
via parallel transport, allows one to map vectors from a tangent space at a point in 
a differentiable curve into the tangent space at another point of this curve. To fix 
notation, consider a differentiable curve a : / —?> M defined in an interval / C R. 
For s,t G /, the parallel transport along a from a(s) to a{t) induced by V is the 
linear isometry denoted by —>■ T^(j.'^M, such that the covariant 

derivative of t e-)- vanishes. If a is continuous and differentiable by parts, its 

parallel transport is constructed joining the corresponding parallel transports along 
each differentiable segment (see e.g. Kobayashi and Numizu |14jL 


2.1 Parallel transport along a curve with jumps 

Let 7 : I —)• M be a cadlag curve with discontinuities in a countable, discrete 
and closed set D = {ti,t 2 , ■ ■ ■}, possibly finite. Suppose that 7 is differentiable in 
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I\D. Let B = (/?n)neN be a family of differentiable curves /3n '■ [0,1] —>■ M such that, 
for each n G N, /3n(0) = lim 'y{tn) and /3n(l) = 7(tn)- The differentiable curves 

B—ytn — 

in the family B fills the gaps along the trajectory of 7 . Hence, we can define the 
parallel transport along 7 with respect to B. Precisely, fix a subinterval [s,t] C I. 
Take J = D (1 {s,t], that is, J is the set of the times that jumps occur in {s,t]. By 
assumption, J is finite and, abusing notation for sack of simplicity, write J = {ti < 
t2 < ■ ■ ■ < tfe}. 

Now, define the curve (7 V B)s,t '■ [s,t + fe] M which concatenates 7 with 
elements of the family {/3i, i = 1,..., A:} in the following way: 


(7 V B)s,t{u) 


/ 

7(u), for u G [s, ti) 

( 3 i{u - ti), for u G + 1 ] 

7(u - 1 ), for u G [ti + l,t2 + 1 ) 

< P2{u -t2- 1), for u G[t2 + l,t2 + 2] 


Pkiu — — k + 1 ), for u € [tk + k — l,tk + k] 

7 (u — k), for u G [tk + k,t + k]. 

Since the constructed curve (7 V B)s,t is continuous and differentiable by parts, we 
define the parallel transport along 7 with respect to B by: 

Note that the domain of {'y\/B)s,t is artificially extended due to the ‘fictitious’ curves 
in B which fill the gaps of 7 . In next sections, the choice of family B will not be 
arbitrary. It will be established by the deterministic flow generated by the vector 
fields of the differential equation. 


3 Stochastic delay differential equations with jumps 

In this section, we present our model of delay differential equations with jumps 
(DDEJ), including the deterministic and stochastic case (SDDEJ). The solution for 
this equation is constructed by induction on the number of jumps, such that, after 
each jump, we use the theory of differential equations without jumps. We remark that 
the existence and uniqueness of solution in stochastic delay differentiable equations 
(without jumps) is a particular case of the theory of stochastic functional differential 
equations (see, e.g. Leandre and Mohammed m)- We start describing the simpler 
context: 


3.1 Deterministic case 

Initially, we construct the discontinuous (cadlag) integrator St which drives our 
model of DDEJ. Let (tn)nGN be an increasing, discrete and closed sequence in R>o 
which indicates the points of discontinuities of St ■ Let (Jn)nGN be the corresponding 
sequence in R of the increments at the jumps of St- Define the integer function 
which counts the number of jumps up to time t hy Nt = max{n : < t}, with 
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the convention that the maximum of the empty set is zero. Consider the integrator 
S : R>o —>■ R in the following way: 


St 


Nt 

t + Jfc. 


fc =0 


The DDEJ in the manifold M is written as: 

dx{t) = Pl,^,{x)F{x{t-d))dSt (1) 

with initial condition given by a differentiable curve /3o : [—d, 0] —M, where d € 
(0,1] is a fixed time delay and T is a smooth vector field in the manifold M. We 
construct a solution ^{t) of equation ([T]) as follows: 

• Solution before the first jump: 

For t G 7 (t) is the solution of the delay differential equation (without 

jumps) given by: 

f x'{t) = F{x{t - d)) 

= Po{t), yt € [-ci,0]. 

• Solution at the jumps: 

Suppose the solution has been constructed in the interval We define the 

solution at the time tm, corresponding to the m-th jump. Consider the ordinary 
differential equations y'n{t) = JuPiUn), for n G N, n < m. We denote the solution 
flows of these equations by (p^. For each n < m, take Zn = lim 7(5). Now, let 

S—¥tn — 

= {ldn)n<m be the family of curves (considered in Section 2) given by: 

Pnit) = Pi (^n)||-Q 

and define 7(tm) = /3m(l)- 

• Solution in the intervals between jumps 

In this case, define the solution using the parallel transport along 7, with respect 
to the family Bm- So, for t G 7(i) is the solution of the following delay 

differential equation: 


f x'{t) = PZfri^) - d)) 

|x(t) = 7(t), Vt G [-d,trn]- 

Note that, although the initial condition of the equation above may have jumps, 
results of the standard theory of delay differential equations on existence and 
uniqueness still hold. In fact, this initial condition is used only to parallel transport 
the vector field, and this (concatenation of) parallel transport has been defined in 
the previous section. 

Therefore, by induction, the solution is well defined for all t > 0. Uniqueness 
comes from the fact that the solution is unique in each step. 

The fictitious curves (as in Section 2) we have used in this construction were 
established by the deterministic flow of the vector field (without delay). For future 
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reference, we call this family of curves associated to 7 by Bp = {Bn,n E N}. This 
model reflects the main motivation that, in some physical situations, the informations 
arriving at a receptor come from different sources (with corresponding different 
delays): here, continuous informations have a fix delay d, but discontinuities in 
the driver integrator have no delay. 


3.2 Stochastic Case 

Let {B}, ..., BY^) be a Brownian motion in a filtered probability space (hi, F, P) 
and {Nt)t>o be a random counting process that indicates the number of jumps up 
to time t, with the properties that Nq = 0 and Nt is hnite (almost surely) for all 
t > 0. Consider a sequence (Jfc)fegN of random variables in Taking B^ = t, 

the integrator of our model is Lt = {L^, Lj,, T™), given by: 


Nt 

L\ = Bl + Y,Jl 

k=0 

An example of this kind of process is the Levy-jump diffusion (see e.g. Applebaum 
[I]), where Nt is a Poisson process and (Jfc) are i.i.d. random variables. 

Write the stochastic delay differential equation with jumps (SDDEJ) by: 

dxt = P^_^ ^A{xt-d) o dLt, (2) 

where Aq, Ai,..., are smooth vector fields in M, considering initial condition 
/3o : [—d, 0] —>■ M, a differentiable curve in M. 

We define the solution of this equation in an analogous way to the deterministic 
case. So, fixing w E fl, in the intervals between the jumps, the solution is given by 
the corresponding Stratonovich stochastic delay differential equation, that is: 

m 

dxt = Y,P7-dd A\xt-d)odLl 

i=0 

with the appropriate initial condition, as we have done in the previous case. Besides 
that, at the times of jump, the solution hops instantaneously in the direction of the 
solution at time one of the following ODE (without delay): 

' m 

' k=0 

y(0) = lim 7(s), 

S-ttn- 


So, with the same notation as before, we have the following result for SDDEJ: 

Theorem 3.1. There exists a unique solution 7 for the SDDEJ ^ defined in t ^ 
[—d, 00), with initial condition 'y{t) = /3o{t) when t E [—d, 0 ]. 

Proof. The existence follows by an analogous construction we have done for the 
deterministic case. Uniqueness holds since that, in each step of the construction, 
the respective solution is unique, by theory of ODE and standard stochastic delay 
differential equations, as in [T7]. □ 
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In this case, for each w G 11, we have a family of differentiable curves established 
by the deterministic flow of the vector held. Again, we call this random family of 
jumps associated to 7 by = {i3„(a;),n € N}. 

The idea of jumping in the direction of the deterministic how at time one comes 
from Marcus SDE, where the integrator is a semimartingale with jumps (for more 
details see Kurtz, Pardoux e Protter m)- 


4 Geometrical aspects of SDDEJ 

In this section, we show that the parallel transport, i.e., the horizontal lift of a 
solution of an SDDEJ in a manifold (M, V) can be described as an SDDEJ in the 
linear frame bundle BM, with respect to a horizontal connection in BM (described 
below). The equation for this horizontal lift corresponds to an extension of the 
results on stochastic geometry started with Ito m and Dynkin [8] to our model of 
SDDEJ (see also [S]). 

4.1 Horizontal lifts to the frame bundle 

For reader’s convenience we recall briehy some geometrical facts about the frame 
bundle of a manifold (for more details, see e.g., among many others, Elworthy [9], 
Kobayashi and Nomizu M)- Let M be a differentiable manifold, with dimension n. 
The frame bundle BM of M is the set of all linear isomorphisms p : 12" TxM for 
X G M. The projection tt : BM —> M maps p to the corresponding x G M. BM is a 
principal bundle over M, with right action of the Lie group GL{n, R), given by the 
composition with the linear isomorphisms. 

Given p in the manifold BM, each tangent space TpBM can be decomposed as 
a direct sum of the vertical and a horizontal subspace, TpBM = VpBM © HpBM. 
The vertical subspace is determined by VpBM = Ker{7r^{p)), where vr* denotes the 
derivative of the projection vr. The horizontal subspace HpBM is established by the 
connection V in M, namely it is generated by the derivative of parallel frames along 
curves in M passing at Tr{p). In this context, one can consider the horizontal lift of a 
vector V G T^M at p G 7 r“^(x) as the unique tangent vector G HpBM such that 

Tr*{p){v^) = V. 

We say that a differentiable curve a : I ^ BM is horizontal when its derivative 
belongs to H^^^f'^BM for all t G /. In fact, given a differentiable curve /? : [0, T) —>■ M, 
and p G 7 r“^(/ 3 ( 0 )), there exists a unique horizontal curve (3^ : [0,T) —)• BM, with 
the property that (t)) = /3(t) for all t in the domain. The curve (3^ is called 
the horizontal lift of the curve /3 (see e.g. [14|L 

Putting together the technique of Section 2 and the horizontal lift described 
above, one can define the horizontal lift of a curve with jumps in M. Let 7 : 
[0, 00) —>■ M be a cadlag curve, D = {ti, • • •} the countable, closed and discrete 
set of points of discontinuity, and B = (/?n)neN be a family of differentiable curves 
I3n '■ [0,1] —>■ M, such that, for all n G N, /3n(0) = lim "/{tn) and /3n(l) = 

S^tn — 

Fix p G 7 r“^( 0 ). Under these conditions, we define the horizontal lift of 7 in p with 
respect to the family B by the cadlag curve in BM: 

7 f^(t) :=P^f{j)op. 
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Each element A in the Lie algebra Ql{n, M) of the Lie group GL{n, R) determines 
a vertical vector field in BM given by, at a point p G BM, 

^*{p) = l*=o- 

The map C//(n,M) i—)• VpBM is surjective. In order to define a SDDEJ in BM, one 
needs a connection in this manifold as well. There are many ways of extending 
a connection V of M to RM. In this section we are interested in the so called 
horizontal lift which is dehned (for a torsion free connection V, see e.g. Cordero 
et al. [ 6 l Chap. 6 ]) as the unique connection in BM which satishes: 


' V%B* = {ABf 
= 0 
= 0 

^ = {VxY)^ 


(3) 


This extension has the property that parallel transport commutes with the horizontal 
lift, that is, if a is a curve in BM, for any v G T^aa{o)Y[^ it holds that PQ^ («)(0 = 
(P(^j( 7 r o a){v))^, see O Lemma 2.1]. 


4.2 Main results 

As in the previous section, initially we deal with deterministic systems. Catuogno 
and Ruffino [31 Prop. 2.1] used the extended connection to prove that, if a 
curve 7 is the solution of a deterministic DDE in M, then its horizontal lift is the 
solution of the corresponding DDE in BM. Now, we extend this result to DDEJ. Let 
F : Af ^ TAf be a smooth vector held and /3o : [—4,0] —Af a differentiable curve 
(initial condition). Let 7 be the solution of the following delay differential equation 
with jumps in M: 


j dx{t) = F{x{t - d)) dSt 

l^x(t) = /3o(t) for t G [—4,0]. 

This solution induces canonically the family Bp oi curves along the deterministic 
how at the jumps, as dehned in Section 3.1. Fix p G 7r“^(7(—4)). Consider : 
[—4,0] BM the horizontal lift of /3o at the point p, and let F^ be the horizontal 
lift of the vector held F. Under these conditions: 


H 13 

Theorem 4.1. The horizontal lift 'fp ’ ^ is the solution of the following DDEJ in 
BM with respect to the connection : 

Idu(t) = Pll,{u) F^{u{t - 4)) dSt 
|u(t) = /3^(t) /or fG [-4,0]. 

Proof. Let (tn)nGN be the increasing sequence of discontinuities of the integrator St. 
Let u : [0, 00 ) —>■ BM be the solution of equation Q, whose existence and uniqueness 
is guaranteed in Section 3.1. We show, by induction on the number of jumps, that 
u{t) = 7 ^’^^(t) for all t > 0. 
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For t € the solution u{t) is given by a delay differential equation, so we 

can apply the result without jumps as in O Prop. 2.1] to obtain the equality. In 
particular, define Pi := lim u{s) = lim 7 ^’^^(s). 

s — s—yti 

At time ti, when the first jump occurs, take the curve /3i G Bp (we recall that 
this is the solution of the ODE y'{t) = JiF{y), with initial condition y(0) = vr(pi)). 
Consider its horizontal lift [3^ at pi. As is the solution of z'{t) = JiF^{z), 
with initial condition 2;(0) = pi, we have that u{ti) = 7p ’ ^(ti). Now, arguing by 
induction, suppose that u{t) = 7p ’ ^{t) for all t G [—d,tm] for m > 1. We claim 
that this equality also holds in the interval In fact, let k be the number 

of jumps that occur in the interval {tm — d,tm)- 

Firstly, for the case {tm+i — tm) > d, we have: 

tm < tm-k + d < . . . < tm-l + d < tm + d < tm+1, 

and it is enough to analyse each of these {k+2) subintervals. For the first subinterval, 
that is, {tm, tm-k + d), we consider a DDE with delay {d + k + 1), where the initial 
condition is the concatenation of the following curves: 

• ut, in the interval [tm — d,tm-k)', 

• Pm-k^ ™ interval [0,1]; 

• Ut, in the interval [tm-k,tm-k+i)', 

• Pm-k+ii ™ the interval [0,1]; 


• Ut, in the interval [tm-i,tm)', 

• Pm, in the interval [0,1]. 

Therefore, for t G {tm,tm-k + d), we have a delay differential equation, with a 
fictitious bigger delay, but without jumps. Applying the result in O Prop. 2.1] the 
equality in this subinterval holds. For the n-th subinterval, with re < {k+2), consider 
a DDE (without jumps), where the initial condition concatenates appropriately ut 
and , with i G {rre — k + n — 1,..., rre} and use [3l Prop. 2.1]. And for the last 
subinterval, that is, for t G {tm + d,tm+i), take the DDE (without jumps) with delay 
d and initial condition given by ut, and again we have the equality. 

Secondly, for the case {tm+i — tm) < d the argument is essentially the same, 
just we have to consider a smaller number of subintervals. Therefore, we have that 
'^(t) = 7p’^^{t) for t G {tm,tm+l)- In particular, define: 

Pm+i ■= lim u{s) = liin 7^’^^(s). 

Finally, consider Pm+i, the horizontal lift of the differentiable curve Pm+i at the 

H 13 

point Pm+i- The equality at tm+i also holds: u{tm+i) = 7p ’ ^{tm+i)- The proof is 
complete. 


□ 


The same result holds in the stochastic case thanks to the transfer principle (see, 
e.g. Emery m), in the sense that; with the same notation as before, let a be the 
solution of the SDDEJ 


I dx{t) = P^_d/x) F{x{t - d))odLt 
=/3o(t) forte[—d,0]. 

This solution induces the random family of curves along the deterministic 

flow at the jumps, as defined in Section 3.2. Hence, we have the following: 

Theorem 4.2. The horizontal lift of a (solution of the SDDEJ in M) is 

the solution of the following SDDEJ in BM, with respeet to the conneetion : 

{du{t) = - d)) o dLt 

/ortG[-d,0]. 

Proof. In the proof of the deterministic case (Theorem 4.1), in each step, the 
solution is given by a standard delay differential equation (without jumps). The 
corresponding lift is the solution of the lifted equation in BM. So, applying 
successively the transfer principle: before the first jump, at the jumps and in the 
intervals between jumps, we have the result (cf. O Thm. 2.2]). 

□ 

The results above do not exhaust the subject; in fact, it opens the possibility of 
exploring the geometry for SDDEJ and corresponding applications: say, holonomies, 
invariant measures, stability (Lyapunov exponents), rotation numbers, etc. 
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